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Polarized parton distribution functions are determined by using asymmetry A\ 
data from longitudinally polarized deep inelastic scattering experiments. From our 
X 2 analysis, polarized n-valence, d-valence, antiquark, and gluon distributions are 
obtained. We propose one set of leading-order distributions and two sets of next- 
to-leading-order ones as the longitudinally-polarized parton distribution functions. 



1 Introduction 

After the EMC finding of a proton-spin issue, many polarized deep-inelastic 
scattering (DIS) experiments have been done on spin structure of the nu- 
cleon. From these experimental data and theoretical studies, we think that 
the nucleon spin is carried not only quarks but also gluons and their angu- 
lar momenta. However, we do not have a clear idea even on the antiquark 
and gluon contributions, which are difficult to be determined by the present 
lepton-nuclcon DIS data. The situation should become clearer in the near 
future because RHIC-Spin experiments will provide valuable information on 
these distributions. 

We tried to determine the polarized parton distribution functions (PDFs) 
by using existing spin asymmetry A\ data for understanding the present situ- 
ation and for suggesting the importance of future experimental studies. The 
following discussions are based on the work in Ref. 1 with the members of 
the Asymmetry Analysis Collaboration (AAC). In Sec. |2|, we explain how to 
calculate A\ in terms of the unpolarized and polarized parton distributions. 
Then, the actual parametrization and x 2 -analysis method are discussed in Sec. 
|j| Our results are shown in Sec. |I]and conclusions are given in Sec. |. 



2 Parton model analysis of polarized DIS data 



There are many measurements of the spin asymmetry A\ for the proton, neu- 
tron, and deuteron. To use these experimental data in our analysis, we express 
Ai as 

9l (x,Q 2 ) 2x[l + R(x,Q 2 )} 



Ai(x, 



Fi(x,Q 2 ) 



F 2 (x,Q 2 



-9i{x,Q 



(1) 
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where F\ and F2 are unpolarized structure functions. The function R(x,Q 2 ) 
is given by R{x, Q 2 ) = <jl/&Ti where <tl and <7y are absorption cross sections 
of longitudinal and transverse photons, and it is determined experimentally in 
reasonably wide Q 2 and x ranges in the SLAC experiment. □ The polarized 
structure function g\{x, Q 2 ) is expressed as 

gi(x,Q 2 ) = ^^eifAC q (x,a 3 )®[Aq i {x,Q 2 )+Aq i (x,Q 2 )]+AC g {x,a s )®Ag(x,Q 2 ) 

i=l 

(2) 

where ej is the electric charge of a quark, and the convolution (g> is defined by 
f(x) ® g(x) = J* &f (^j g(y). The distribution A% = qj - q\ represents the 
difference between the number densities of quark with hclicity parallel to that 
of parent nucleon and with hclicity anti-parallel. The definitions of Aq~i and 
Ag are the same. AC q and AC g are the coefficient functions. In discussing 
unpolarized reactions, the structure function F2 is usually used rather than 
Fi, and F2 can be written in terms of unpolarized PDFs, qi, q~i, and g, with 
coefficient functions in the similar way to g\ . In the next-to- leading-order 
(NLO) analysis, we choose the modified minimal subtraction (MS) scheme. 

We provide the polarized parton distributions at Q 2 = lGeV 2 (= Qq). 
Then, the distributions are evolved from Q 2 to experimental Q 2 points by 
DGLAP equations. In our numerical analysis, we use a modified version of the 
program in Ref. 3, where the evolution equations are solved by a brute- force 
method. 



3 Parametrization of polarized parton distributions 

Now, we explain how the polarized parton distributions are parametrized. The 
unpolarized PDFs fi(x,Qg) and polarized PDFs Afi(x, Qq) are given at the 
initial scale Qq. Here, the subscript i represents quark flavors and gluon. 
In our analysis, we require the positivity condition of the PDFs in order to 
constrain the forms of the polarized PDFs. Therefore, it is convenient to take 
the following functional form of the polarized PDFs at Qq: 

Afi(x,Q 2 Q ) =A lX a > (l+7 4 x A ')/,(x,Q2). (3) 

The positivity condition is originated in a probabilistic interpretation of the 
parton densities: the polarized PDFs should satisfy | Afi(x, Q 2 ) | < fi(x,Q 2 ). 
In our analysis, we simply require that this condition should be satisfied not 
only in the leading-order (LO) and but also in the NLO at Qq. Thus, we have 
four parameters (A4, a.i, 7^ and A^) for each i. 
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In addition to the positivity condition, wc assume the SU(3) flavor sym- 
metry for the sea-quark distributions at Q\ to reduce the number of free pa- 
rameters. Then, the first moments of Au v (x) and Ad v (x), which are written 
as T] Uv and r\& v , can be described in terms of axial charges for octet baryon, F 
and D, measured in hyperon and neutron /3-decays. They are determined as 
F = 0.463 ± 0.008 and D = 0.804 ± 0.008, which lead to r) Uv = 0.926 ± 0.014 
and r/d v = —0.341 ± 0.018. In this way, we fix these first moments, so that 
two parameters A Uv and Ad v are determined by these first moments and other 
parameter values. Thus, the remaining job is to determine the values of the fol- 
lowing 14 parameters, Aq, A g , a>i, 7j, Xi (i = u v , d v , q, g), by the x 2 analysis 
of the polarized DIS experimental data. 

We determine the values of the 14 parameters by fitting the Ai(x, Q 2 ) data 
for the proton from E130, E143, EMC, SMC, and HERMES, the neutron from 
E142, E154, and HERMES, and the deuteron from E143, E155, and SMC. 
We also use LO and NLO GRV parametrizations for the unpolarized PDFsa 
and the SLAC measurement of R(x,Q 2 ).a Then, the best parametrization is 
obtained by minimizing 



(Af ata (x, Q 2 ) - Af^(x, Q 2 )) 2 
(AAf ata (a;,Q 2 )) 2 



(4) 



where AA^ ata represents the error on the experimental data including both 
systematic and statistical errors. In evolving the distribution functions with 
Q 2 , we neglect charm-quark contributions to Ai(x,Q 2 ) and take the flavor 
number Nf = 3, because the contribution is very small in a few Q 2 region 
where most experimental data exist. To be consistent with the unpolarized, 
we use the same scale parameters as the GRV,Q Aqq D = 204 MeV at LO and 



A, 



(3) 
QCD 



299 MeV at NLO. 



4 Results 

We got minimum \ 2 - x 2 /d.o.f=322. 6/360 
for the LO and x 2 /d.o.f=300.4/360 , 
for the NLO. The difference between 
the LO and NLO x 2 values is about 
7%, which indicates the importance 
the NLO analysis. 

We show the Q 2 dependence of 
spin asymmetry A\ for the proton in 
Fig. 0. This figure indicates that the 
NLO effects become larger in the small 
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Figure 1: Theoretical of spin asymmetries A\ 
for the proton are compared with the exper- 
imental data at x « 0.117. (from Ref. 1) 
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Q 2 region and that there is strong Q 2 dependence especially in the small Q 2 
region. It is not right to assume Q 2 independence of the spin asymmetry Ai 
in obtaining g\, so that, we have to be careful using PQCD in this region. 

Next, we show the polarized par- 
ton distributions for the NLO in Fig. 
||. As the figure indicates, we obtain 
negative polarization for Aq(x), and 
large positive polarization for Ag(x). 
The first moment for Au v (x) is fixed 
at the positive value and the one 
for Ad v (x) is fixed at the negative 
value, so that the obtained distribu- 
tions Au v (x) and Ad v (x) become pos- 
itive and negative, respectively. Simi- 
lar results are obtained in the LO dis- 




Figure 2: Obtained NLO polarized parton 
distributions at Q 2 = 1 GeV 2 . (from Ref. 1) 



tributions.Lll There are slight differences between the LO and NLO distributions 
Ad v {x). However, the differences are large between the LO and NLO gluon 
distributions in the wide x region. It is caused by the gluon contribution 
through the coefficient function. We calculate the quark spin content by using 
the obtained LO and NLO distributions. It is given by AS = r] Uqj + rjd v + 6 r]q, 
where rjq is the first moment of Aq(x). Because rj Uv and are fixed, only r\q 
affects the spin content in the different analyses. The LO and NLO moments 
are r]q — —0.064 and —I). 089, so that the spin content becomes AS = 0.201 
and 0.051, respectively. El The NLO spin content (AS = 0.051) is significantly 
smaller than other analysis results. For example, the recent SMC and Leader- 
Sidrov-Stamenov (LSS) parametrizationstm obtained AS =0.19 and 0.28 at 
Q 2 = l GeV 2 . In order to investigate the reason for the small AS in our analysis, 
we show each antiquark distribution in Fig. |^. 

The NLO antiquark distributions of 0J " 
the SMC, LSS, and AAC analyses are °- 02 ; 
calculated at Q 2 =1 GeV 2 . Because 0M 
the antiquark distribution is not directly 
given in the SMC analysis, we may call 
it as a transformed SMC ("SMC") dis- 
tribution. It is calculated by transform- 
ing the published distributions by the 
SMC. All the distributions agree in prin- 
ciple in this region (0.01 < x < 0.1) 
where accurate experimental data exist 
and the antiquark distribution plays an 



— AAC 
LSS 

— "SMC" 



Figure 3: Antiquark distributions of recent 
parametrizations "SMC" , LSS and AAC at 
Q 2 =l GeV 2 . (from Ref. 1) 
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important role. On the other hand, it is clear that our distribution does not fall 
off rapidly as x — » in comparison with the others. This is the reason why our 
NLO spin content is significantly smaller. In fact, we obtained the parameter 
a for the antiquark distribution as aq(NLO) = 0.32 ± 0.22, which controlled 
the small- a; behavior of Aq(x). However, the large error of the parameter aq 
suggests that the small- a; part of Aq(x) cannot be fixed by the existing data. 
Actually, there is no data in the small- a; region (x < 0.04). Therefore, we had 
better consider to constrain the parameter aq by theoretical ideas. We discuss 
such possibilities by using the Regge theory and the perturbative QCD. 

According to the Regge model, the small- a; behavior of g\ is suggested as 
gi(x) ~ x- a with a = -0.5 ~ 0.0.0 Therefore, we expect Aq(x) ~ x -°~ - 5 
as x — > 0. Because the parametrized function is given by Aq(x)/q(x), we 
should find out the small- a: behavior of the unpolarized distribution. The 



GRV distribution has the property x q 



at Q 2 —l GeV 2 according to 



our numerical analysis, the Regge prediction becomes a^ e " e = 1.1 ~ 1.6, if 
the theory is applied at Q 2 — 1 GeV 2 . The perturbative QCD could also 
suggest the small- a; behavior. If we can assume that the singlet-quark and 
gluon distributions are constants as x — > at certainX) 2 (= Q\), their singular 
behavior is predicted from the evolution equations. The singlet distribution 
behaves like AS(x) ~ x a as a = —0.12 ~ —0.09, if we choose the evolution 
range from Q\ = 0.3 ~ 0.5 GeV 2 to Q 2 = 1 GeV 2 . Therefore, the perturbative 



vQCD 



1.0. 




< O.li* 06=05 



0.00001 0.0001 



QCD with the assumption of the above Qf range suggests a| 

The Regge theory and perturbative 
QCD suggest the range aq = 1.0 ~ 1.6, 
so that we try the NLO x 2 analyses by 
fixing the parameter at aq =0.5, 1.0, 
and 1.6. The last two values are in 
the theoretical prediction range, and the 
first one is simply taken as a slightly sin- 
gular distribution. The obtained mini- 
mum x 2 values are larger than the NLO 
fit (x 2 =300.4) by 0.1, 1.8, and 7.7%, and 
the first moments are A£=0.123, 0.241 
and 0.276 for aq =0.5, 1.0, and 1.6, respectively. The small- a; falloff for larger 
aq changes the rjq and AE significantly. We show the spin content in the 
region between x m i n and 1 by calculating AS = AT,(x)dx in Fig. |j. 

Because the LSS and SMC distributions are less singular functions of x, their 
spin contents saturate even at x m i n = 10~ 4 although the AS of our NLO re- 
sult with free aq still decreases in this region. If the parameter aq is taken in 
the perturbative QCD and Regge theory prediction range, the calculated spin 



Pi gure 4: ^rain dependence of AS at Q 
lGeV 2 . (from Ref. 1) 
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content is within the usually quoted values AE = 0.1 ~ 0.3. In this sense, our 
results are not inconsistent with the previous analyses. Our results indicate 
that the spin content cannot be determined uniquely, because the accurate 
experimental data are not available in small x region. The obtained x 2 value 
suggests that the ag=1.0 solution could be also taken as one of the good fits. 
The ctq=0.5 distributions are almost the same as the ones in the free-ckq NLO 
analysis, so that it is redundant to propose it as one of our good fits. 

From our analyses, we propose the LO distributions, the NLO ones with 
free Uq (NLO-1), and those with fixed a.q=l.O (NLO-2) as the longitudinally- 
polarized parton distributions of the A AC analyses. Useful functional forms 
are given at Q 2 =l GeV 2 in Appendix B of Ref. 1 for practical applications. 

5 Conclusions 

From the LO and NLO % 2 analyses, we obtained good fits to the experimen- 
tal data. Because the NLO x 2 is significantly smaller than that of LO, the 
NLO analysis should be necessarily used in the parametrization studies. It is 
particularly important for extracting information on Ag. However, the polar- 
ized antiquark and gluon distribution cannot be uniquely determined by the 
present DIS data. We provide the optimum LO and NLO distributions at 
Q 2 = 1 GeV 2 from our numerical analyses. 
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